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We study the electronic instabilities in a 2D Hubbard model where one of the dimensions has a 
finite width, so that it can be considered as a large array of coupled chains. The finite transverse 
J"""""- ' size of the system gives rise to a discrete string of Fermi points, with respective electron fields that, 

due to their mutual interaction, acquire anomalous scaling dimensions depending on the point of the 
string. Using bosonization methods, we show that the anomalous scaling dimensions vanish when the 
number of coupled chains goes to infinity, implying the Fermi liquid behavior of a 2D system in that 
limit. However, when the Fermi level is at the Van Hove singularity arising from the saddle points of 
the 2D dispersion, backscattering and Cooper-pair scattering lead to the breakdown of the metallic 
behavior at low energies. These interactions are taken into account through their renormalization 
group scaling, studying in turn their influence on the nonperturbative bosonization of the model. 
We show that, at a certain low-energy scale, the anomalous electron dimension diverges at the Fermi 
points closer to the saddle points of the 2D dispersion. The d-wave superconducting correlations 
become also large at low energies, but their growth is cut off as the suppression of fermion excitations 
takes place first, extending progressively along the Fermi points towards the diagonals of the 2D 
I ■ Brillouin zone. We stress that this effect arises from the vanishing of the charge stiffness at the 

Fermi points, characterizing a critical behavior that is well captured within our nonperturbative 
CZ3 , approach. 
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Q ; I. INTRODUCTION 

Copper-oxide superconductors have shown since their discovery quite unusual electronic properties, regarding the 
particular features of the superconducting state as well as the behavior of the normal state. One of the most remarkable 
experimental observations has been the existence of the so-called pseudogap, at temperatures which are above the 
superconducting dome in the phase diagram^. The unconventional character of this pseudogap lies in that it opens up 
quite anisotropically along the Fermi line of the electron system. At some doping-dependent temperature T*, which 
is in general well above the transition temperature of the superconducting phase, electron quasiparticles become 
suppressed in the region of the Fermi line which is closer to the momenta (tt : 0) and (0, n). This is particularly 
manifest in the ARPES experiments^, as the quasiparticle peaks get a dispersion that does not cross the Fermi level 
at any point in momentum space. As the temperature is lowered from T* , this partial destruction of the Fermi line 
takes place progressively over wider regions, approaching the diagonals of the Brillouin zone^. 

There have been several attempts to provide a theoretical explanation for the pseudogap in the cuprates. Some of 
the proposals have considered the possible coupling of the electron quasiparticles to collective boson excitations^. 
It has been also investigated the possibility that the electron system may develop an ordered phase, which would 
induce the appearance of a gap in the single-particle spectrum^. Other explanations are based on the strong electron 
correlations and, in particular, on the description of the cuprates in terms of the Hubbard model. Thus, it has been 
proposed that theories based on a resonating-valence-bond state 8 or on the existence of states with time-reversal 
symmetry-breaking^ may give rise to features consistent with the pseudogap phase. 

The main challenge that face all these proposals is actually the difficulty to deal with a regime of strong electron 
correlations. In these conditions, one has usually to resort to mean-field approximations, that may give a qualitative 
idea of the possible phases of the electron system. It would be also desirable to have an alternative framework in 
which the nature of the pseudogap feature could be clarified. In this regard, the motivation of the present paper is to 
investigate the Hubbard model by means of a nonperturbative approach, trying to capture the effects of the strong 
electron correlations. For that purpose, we will consider a 2D model where one of the dimensions has a finite width, 
so that it may be considered as formed by the coupling of a large number of Hubbard chains. A general analysis of 
the iV-chain Hubbard model in weak coupling has been accomplished in Ref. [l(J. We will use bosonization techniques 
to determine the nonperturbative effects of the interactions with small momentum-transfer on the properties of the 



(N 



O 



2 



fermion excitations. In this framework, the Fermi liquid properties of a 2D system are recovered when the number of 
coupled chains goes to infinity, as the anomalous scaling dimensions of the electron fields vanish in that limit. The 
description has to be completed anyhow by considering the interactions whose momentum-transfer is not small. These 
can be taken into account through their renormalization group scaling, studying their influence on the nonpcrturbative 
bosonization of the model and the instabilities they may trigger at low energies. 

The most interesting instance corresponds to the case where the Fermi level is near the Van Hove singularity arising 
from the saddle points in the 2D dispersion at (tt, 0) and (0,7r). When the 2D model is approached by a large array 
of coupled chains, the fermion excitations are attached to a large but finite number of Fermi points, with respective 
electron fields that are renormalized by the interactions and characterized by their anomalous scaling dimensions. We 
will see that, for the model with on-site repulsion, the angle- resolved quasiparticle weight vanishes below a certain 
low-energy scale, starting from the Fermi points which are closer to (tt,0) and (0,tt). The <i-wave superconducting 
correlations are also dominant at low energies, but their growth is cut off as the suppression of fermion excitations 
takes place first, extending progressively along the Fermi points towards the diagonal of the Brillouin zone. 

The results that we present are in part reminiscent of those obtained in 2D electron systems when the Fermi line is 
placed close to the saddle points at (tt, 0) and (0, Tr^dLl^JA^J^IJ^J^^l^^l. It has been shown that the presence 
of a Van Hove singularity induces a strong renormalization of the quasiparticle properties, with a clear suppression 
of the quasiparticle weight near the saddle point s 15 ! 23 . Some studies have also remarked that the divergent flow of 
the couplings observed in renormalization group analyses could lead to the vanishing of both the compressibility and 
the magnetic susceptibility 1 ^ 19 ' 20 . In our approach, we are able to identify the onset of the destruction of the Fermi 
line at an energy scale where the electron correlations have not entered yet the divergent regime. This is possible as 
long as the Fermi surface instability appears as a critical point in the nonperturbative quasiparticle renormalization 
from the bosonization approach. The present combination of bosonization and renormalization group methods may 
be seen as an alternative to conventional analyses of 2D systems near the Van Hove filling, offering a sensible 2D limit 
as the number of coupled chains is increased. 

II. QUASIPARTICLE PROPERTIES IN THE LIMIT OF LARGE NUMBER OF CHAINS 

We will take as our starting point a 2D model with nearest-neighbor and next-to-nearest-neighbor hopping, with a 
finite size in one of the spatial dimensions. For the time being, we introduce a generic spin-independent interaction 
with potential V. The hamiltonian of the system is then given by 

N oo 

H = (-*4r(n)iMn + l)-tl&(»)l^^ 

j—l n— — oo 
N oo 

+ Y Y ^jffM^M V#(m,n) CW^'W (!) 
j,l=l m,n— — oo 

where ipj a ( n ) anc ^ 1 Pj^( n ) are electron creation and annihilation operators at site n of the j-th chain, with spin 
projection given by the a index. 

For the description of quasiparticle properties, it is convenient to make the passage to electron operators in mo- 
mentum space, taking the Fourier transform in the site and chain variables: 

^An) = 4= E f d P e-^e-^ a (k a ,p) (2) 

V ^ a J-7V 

We will apply periodic boundary conditions to the system of N chains. Whenever it is convenient to exclude the 
singular points (tt, 0) and (0, 7r) from the spectrum, we will make however a twist by e l7r in the electron field after 
completing a whole period along the array. With these twisted boundary conditions, the transverse momentum k a 
will be quantized according to the rule k a — 7r(2a + 1)/N, a — —N/2, —N/2 + 1, . . . N/2 — 1. From the hamiltonian 
([T]), we get the dispersion relation for the different transverse momenta 

£ a(p) = —%t cos(p) — 2t cos(fc a ) — 4i' cos(p) cos(fc a ) (3) 

Furthermore, for the sake of describing low-energy properties, we may concentrate on the modes that are within a 
cutoff energy E c about the Fermi energy ef, where the dispersion becomes approximately linear about the Fermi 
points. This is illustrated in Fig. [TJ For each transverse momentum k a , the low-energy excitations can be encoded 



into two fields ij) aa {p) and ipatrip) with opposite chirality (i.e. right- or left-moving character), corresponding to the 
two different Fermi points at p a and —p a : 

dpipl(k a ,p) {s a (p) - e F )i> a (k ai p) 

dpv a {i>\ a (p)(p ~ p a )i> a<y {p) + 1pljp)(-P - Pa)^ha(p)j (4) 

v a being the Fermi velocity at p a and — p a . 

(0,71) 
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FIG. 1: Schematic representation of the Fermi points for a 2D model made of a finite number of coupled chains, when the Fermi 
level is at the Van Hove singularity. The horizontal lines correspond to quantized transverse momenta for twisted boundary 
conditions, k a = n(2a + 1)/N, with a = —N/2, —N/2 + 1, . . . N/2 — 1. Energy contour lines are also represented for the cutoff 
at E c and — E c about the Fermi level. 

From the detailed discussion in Ref. [l(| it follows that E c can be taken suitably small so that the possible scattering 
processes are severely limited by momentum conservation. The price that one has to pay for the reduction of the 
cutoff is that the original interactions are corrected by terms which are of the order of ~ (V/2irt) 2 log(kF / E c ), kp 
being the scale of the size of the Fermi line. In our model, we will start with an on-site interaction U such that 
U/2irt < 1, in such a way that the perturbative corrections are under control down to a conveniently small cutoff E c . 

Confining our discussion to the low-energy modes in the interval between sf — E c and ep + E c , there is a reduced 
number of ways in which the momentum can be conserved in the scattering processes. These can be classified into 
four different interaction channels, as represented in Fig. [21 We may have processes with suitably small momentum 
transfer between interacting particles near Fermi points a and b (to which we will assign couplings /i^ ) and processes 
where the particles exchange their positions in the neighborhood of the two Fermi points (which we will label with 
couplings fab )- The other type of interactions that exist generically at low energies corresponds to processes where 
the total momentum of the incoming particles is almost vanishing. In this case, we can discern a channel where the 
particles interact keeping their respective chiralities (with couplings that we will denote by c£j ) and another where 

they exchange their chiralities (requiring then different couplings }^r^- The hamiltonian involving the low-energy 
excitations becomes then 

H = Yl / d P Va { 1 PL(P)(P - Pa)^Pacr(p) + 

+ J2 I d P d P' d 1 [ {4a (P + <Z)^W (p) /$ V»L ^ ~ OWl* ( P ') 
<a,b> J 

+4Ap + <?)Vw(p) /# 4Ap - q)A*'(p') 

+4 la (P + ?)Vw (p) f { ~^ 4>t a > (p' ~ ip') + a a, b <-> bj 

+4 b JP + l)^a*(p) ci+ } ^l h AP' - 3)^-00* (pO 

W- h „{p + ?)VW(P) 476 4a' (P' - «)V-a^ (P )} (5) 



where the indices in the sum over pairs < a,b > run over just one of the chiralities. 
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To clarify our notation, we recall that the index a corresponds to a Fermi point with the same transverse momentum 
but opposite chirality to that represented by the index a. Thus, the difference between the interactions ff^~ and 

is that the first involves modes which have the same chirality, while / "t couples an electron current made of right- 

a,b 

moving fields to another current of left-moving fields. On the other hand, the difference between the interactions 
represented by / ^ and - is that, in the latter, the incoming particles exchange their chirality in the scattering 

a.b a,b 

process. We finally remark that, in our notation, the index —a represents a Fermi point with a momentum completely 
opposite to that of Fermi point a (so that the modes around Fermi point —a have opposite chirality to those around 
Fermi point a). In this way, the incoming particles are bound to have very small total momentum in the interactions 

represented by and cj^— . 




(a) (b) 





FIG. 2: Interaction processes that survive at low energies within the cutoff E c about the Fermi level. The solid and dashed 
lines correspond to the two different chiralities (right- and left-moving character) of the electron fields. 

An important difference between the interactions of particles with like chirality (ff^, /_- , f^b^ /- ) an< ^ those 



corresponding to / % , f y J , c { V , c} n J is that the strength of the latter depends on the energy scale. This may induce 

a, b a,b ■ ■ 

instabilities that appear when the system is observed at scales very close to the Fermi level. On the other hand, the 
significance of the interactions with small momentum transfer , /_- , / - , fjj^) is that they dictate the class 

of electron liquid describing the normal state (metallic) properties of the coupled chains. This corresponds to the 
so-called Luttinger liquid for any finite number of chains^. We discuss next the way in which the limit of a large 
number of chains may change the quasiparticle properties, depending on the given Fermi point. 

The Luttinger liquid properties are best described in terms of operators given by the density fluctuations at small 
momentum: 



\/\q\ J 



q>0 



4(9) 



dp i>\a (P + 1)^ (P) 1 < 



(6) 
(7) 



It is well-known that the kinetic part of the hamiltonian can be written in terms of the operators ([6]) and J7| 26 . 
For the time being, we will restrict the discussion to the interactions preserving the chirality of the particles. Then, 
we end up with the hamiltonian governing the Luttinger liquid dynamics 



H 



LL 



9>0 



dq \q\ B\{q)B a {q) + / dq \q\ B^B^q) 



q<0 
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+ E / *|g|(/iJ(^(?)56(?) + Ba(g)B 6 t (g)) + /^ ) (Bt(g)^(-g) + fl a (g)^(-g))) 

<a,b> ^9 >0 

+ E / d g|gl(4V(4(g)^(«) + ^K^)) + f^(Bl(q)Bl(-q) + Mq)B b (-q))) (8) 

<a : 6>- 7 '?< V ' 

The hamiltonian can be diagonalized by making the change of variables 

Blid) = E(^4(?)+^(-<z)) (9) 
& 

Ba(-q) = Y,( S ab B b(-l)+ t abBt( q )) (10) 



The change of basis defined by the matrix elements s a b, t a ^ is not unitary, since the requirement of canonical commu- 
tation relations for the new operators imply in particular that 

E(M 2 ~l^l 2 ) = i (ii) 

b 

For the sake of carrying out the diagonalization of Hll with great computational efficiency, one can introduce however 
the following trick. By redefining all the boson operators with left-handed character, for instance, 

Ml)=iMq) , Bl(q)=iBl(q) (12) 
we change the canonical commutation relations to 

%(q),B}_(q')] =-6(q-q') (13) 

Keeping B a (q) = B a ( q ), B\( q ) = B\( q ) for the other chirality, the linear change of variables for the redefined operators 
has matrix elements that correspond to a unitary change of basis. Thus, by performing first the redefinition (|12[) . one 
can apply then standard algorithms that are able to diagonalize very large complex matrices. By reverting at the end 
the transformation (fl"2")) , it is possible then to bring the hamiltonian ([5]) to the form 



H LL =J2va([ dq\q\Bl(q)B a (q)+ [ dq\q\Bl{q)%{q) 

\Jq>0 Jq<0 



(14) 



Our aim is to show the evolution of the quasiparticle properties at the hot spots around the saddle points (7r, 0) and 
(0, 7r). As a most relevant instance, we give the results of performing the diagonalization of the hamiltonian ([8]) for a 
growing number of coupled chains with t' = —0.2t, at the filling that makes the 2D Fermi line to go across the saddle 
points. We have taken an on-site interaction U , that leads to a set of momentum-independent couplings equal to U/N 
in the hamiltonian ([5]). In general we may interpret that the original boson operators are dressed by the interactions. 
The different correlation functions can be obtained by expressing the operators B\{ q ) and B a ( q ) in terms of the set 
of noninteracting operators in (|14[) . in such a way that the interaction effects are encoded in the string {v a } and in 
the parameters of the transformation (jTD|) . 

The most important piece of information comes, after the diagonalization, from the change of basis that brings the 
hamiltonian to the form (|14[) . The original boson operators can be obtained from the free operators that appear in 
(fill through the inverse change of variables 

B l(D = E( S ^(9)+^4(-<?)) (15) 

b 

M-q) - J2{ s ab B b(- ( i) + ^bBl( q )) (16) 

6 

In the absence of interaction, the only nonvanishing coefficients are the diagonal elements s aa — S^a = 1- The 
interactions give rise however to anomalous scaling dimensions 7 a of the electron fields, which can be obtained from 
the spatial decay of the electron propagator at equal time t 

(i>l(x,t)^ aa (y,t)) ~ — i— - (17) 

\ x y\ 
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FIG. 3: Dependence of the anomalous electron dimension j a on the number N of coupled chains in a Hubbard model with 
periodic boundary conditions, when the Fermi point is taken close to the diagonal of the Brillouin zone (full line), with a 
constant transverse momentum k a — 7r/64 (dashed line), and at the nearest location to the saddle point with k a = 2n/N 
(dotted line). 



It is well-known that the fermion fields Vw, V'ao- are recovered by exponentiation of the integral of the corresponding 
chiral part of the uniform density p a (x)^. This is given according to ^ and ([7]) by 

Pa{x) « / dq^Wl ^ (Bt(q) + Bsi-q)) + [ dq^~\ e"^ (4(g) + B a (-q)) (18) 

Jq>0 Jq<0 V 7 

Then, it can be seen that 

1 + 7 a =E(l^! 2 + l^| 2 ) (19) 

b 

We observe that, as soon as the interaction is switched on, we get an anomalous dimension j a ^ 0, reflecting the 
disappearance of the quasiparticle pole in the electron propagator. 

It is very instructive to observe the behavior of the anomalous dimension 7 a at different Fermi points of the model. 
We find that, away from (jr, 0) and (0, 7r), the anomalous dimension -f a decreases in the limit of large number N of 
chains, following with great accuracy a 1/N behavior—. The prefactor for this law depends however on the distance 
to the saddle points, and it grows as one gets closer to (tt, 0) or (0, 7r), as shown in Fig. [3] The limit towards 
these points becomes actually quite singular. If one takes the limit of large N sitting at the Fermi point which is 
away but nearest to the saddle point (so that its Fermi velocity decreases consequently as 1/N) we observe that the 
anomalous dimension also decreases monotonically, but with a behavior that seems to be smoother than any power- 
law, as observed from Fig. [3j The sizes for which we have been able to diagonalize the model (up to w 2000 chains) 
do not allow us to assure whether this limit to the saddle point may produce a vanishing anomalous dimension as 
N — ► oo. Our results make clear anyhow that, for a Fermi point placed at any finite distance from the saddle points, 
the anomalous dimension vanishes in the limit of very large number of chains. This is what allows to take the limit 
N — > co as a way of approaching the properties of an interacting Fermi liquid in two dimensions. We will see actually 
that the presence of backscattering and Cooper-pair scattering modify significantly the present picture, leading to an 
unconventional low-energy behavior near the hot spots for large number of Fermi points. 



III. RENORMALIZATION GROUP APPROACH TO LOW-ENERGY PHASES 



In order to describe the quasiparticle properties at low energies, one has to incorporate the effects of the variation 
of the couplings f r j c i + & i c i b with the energy scale. The discussion of the preceding section is pertinent at the 
energy scale E c where the dispersion becomes approximately linear about the Fermi level. At lower energy scales, 
the effective values of the couplings may change significantly, having a large impact on the renormalization of the 
quasiparticle properties. 

We have to remark at this point that the renormalization effects manifest in general not only in the quasiparticle 
parameters, but also giving rise to changes in the shape of the Fermi surface. It has been proposed that, within 
a renormalization group approach at constant particle number, the Fermi surface has to be understood as a fixed 
point of the renormalization group2£. A practical implementation of this idea for a two-chain model can be found 
in Ref. |2jJ. In the present context, to find such a fixed-point solution is beyond reach for the large string of Fermi 
points, as the Fermi surface is reshaped also by irrelevant interactions whose momenta are not constrained by the 
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cutoff, and that proliferate at large N even more rapidly than the couplings already considered. Being interested in 
the situation where the Fermi level is close to a Van Hove singularity, it will be enough for us to guarantee that the 
Fermi line is not driven away from the saddle points at (tt, 0) and (0, tt) upon switching the interaction. We rely for 
that on the evidence obtained with the renormalization group at fixed chemical potential, which provides a sensible 
approach to systems in contact with a large reservoir of particles. In this context, it has been actually shown that the 
Van Hove singularity acts as an attractor, pinning the Fermi level for a certain window in the choice of the chemical 
potentia l 12 ! 15 . We motivate in this way the suitability of studying the model with the Fermi line near (tt, 0) and 
(0, 7r), as this choice does not imply the fine-tuning of the chemical potential in the system (as fixed for instance by 
the contact with a charge reservoir). 

Below the scale E c , the four-fermion interaction vertices are corrected by diagrams that depend logarithmically on 
the energy scale A (measured from the Fermi energy £p). To deal with this divergence, one can integrate progressively 
the electron modes starting from the limits of the linear branches at £p ± E c 30 . This leads to a dependence of the 
couplings on the scale I = \og(E c /A), governed by scaling equations that read at the one-loop level 

d 4f = -^ b ^ 2 -^ ^ 



Of 



(-) 



a > b - 1 ffA-)\2 , f „{-) ^^2 „(-)„(+) 



dc {+) 



dl TTVab^b' ^^a,-b J a-ba-b> V ' 

81 ^ 2ttv c [Ca > c Cc - fc + ° a > c Cc < b ' + 2nv ab a '» a >- b 1 ' 

C 



dl ^ 2ttv, 



c <^c,b " r C a.c c c.b 



l-cC-M-) LVc'^'^ (23) 

C a,b' l a,-b r, , /, C a,b' l a,-b V ZJ 7 
TTVab 2TTV ab ^— ' ' 

s 

where h^_ b = 2f^_ b — S a bC^ a , and the v a i, are defined in terms of the Fermi velocities as v a b = (v a + Ub)/22i. 

We have solved numerically the renormalization group equations for a Hubbard model with next-to-nearest neighbor 
hopping t' = — 0.2i. Our main interest focuses on the instabilities arising from the hot spots around the momenta (ir, 0) 
and (0, 7r). Thus we have chosen the filling to match the Fermi level with the position of the Van Hove singularity from 
the saddle points. We have considered systems with up to N = 76 and 152 Fermi points, which already demand the 
introduction of « 5900 independent couplings. We have assumed that the system becomes periodic in the transverse 
direction, but using twisted boundary conditions to avoid the presence of Fermi points with vanishing Fermi velocity 
at (tt, 0) and (0, 7r). Anyhow, in this approach we are able to capture the physics of the hot spots by increasing the 
number N, since the lowest Fermi velocity (corresponding to k a = ir/N and k a = tt — tt /N) scales as 2n{t 2 -At' 2 ) 1 / 2 /N . 
It is worth to point out that the model remains exactly C^-symmetric with the mentioned choice of filling level. 

We have started the integration of the equations in the weak coupling regime, taking a value of the on-site Hubbard 
repulsion U — t. We recall that, when making the passage to the electron fields in momentum space, all the couplings 
in the hamiltonian ((5|) get a contribution equal to U/N from the on-site interaction. In any event, the flow is driven 
to a regime where some of the couplings grow large. In particular, at some energy scale A* = E c e~ l , a divergence 
is found in some of the couplings. For N = 76, for instance, we find that I* — 4.99. As observed in Ref. [l0, the 
singularity corresponds to the divergence of interactions involving scattering between the Fermi points with smallest 
velocity, i.e. those nearest to the hot spots. The rest of the couplings remain in the weak coupling regime. 

In order to characterize the instability from the large growth of the couplings, we have also computed the scaling 
of the response functions R^\uj) for singlet (SS) and triplet superconducting (TS) order parameters, as well as for 
charge-density-wave (CDW) and spin-density-wave (SDW) order parameters. The superconducting order parameters, 
for instance, are given by 



J dp(* aT (p)f_ Qi (-p) - tf oi (p)*_ oT (-p)) (24) 



Oi TS) = J dp(*« T (p)*-ai(-p) + *«i(p)*-aT(-p)) (25) 
These operators build response function matrices defined in the space of Fermi points 

f+OO 

4%) - -i / dte^(TO^(t)O^(0)) (26) 
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with (rj) = (SS),(TS). 

Due to the nontrivial structure of the Fermi line, the symmetry of the pair wavefunction is obtained by diagonalizing 

(SS) (TS) 

the matrices R a b and R ab in the low-energy limit, and looking for the eigenvector corresponding to the largest 
eigenvalue. We recall at this point that the derivatives with respect to the frequency of the response functions, 

R = dR^ /d\og(uj), have well-defined scaling properties^. In the case of the superconducting response functions, 
the scaling equations read: 

(55) 

= ^ Eta - y^' + 47c ) ) + 5aM ( :i + /UpS s) m (27) 

(T5 1 ) 

Because of the symmetry of the Fermi line, the eigenvectors of R a fi transform according to the irreducible represen- 
tations of Ci v . 

In the case of density-wave response functions, given a momentum-transfer q across the Fermi line connecting the 
Fermi points (a, b), the only other pair connected by q is made of the inverted points (—6, —a). We may form therefore 
two different operators for the same momentum q 

[ CDW/SDW \ q ) = Jdp(^(p + q )^(p)±^ al {p + q ^- bl (p)) (29) 

(CDW/SDW) {ci) = y'^^t__ T(p + g ) ^_ oT(j3)± ^ (p + ^_ oi(p) ) (30) 

which give rise to two different response functions 

/+oo 
dte- t <TO^ ) (q,t)oW t (q,0)) (31) 
-GO 

/+oo 
^e- t <TOW(q,t)oW t (q,0)) (32) 
-oo 

The response functions i?^(q, uj) and R^ i \{(\ 1 lo) are actually entangled by the interactions, as can be seen from 
the scaling equations 

a -j-(CDW) 

1 '"•" ) - [«& - ^rv^ + (4? - ^sr-wi] (») 

<« - 2e<-> ^"''(q,-) + (& - 2/<-')RS D '"'(q^)] (34) 
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aiog(A) 


7r(u 4 


■ v b ) 



_ b l 1,1 ; w a>fc ■> ab 

(+) — (5TW) , (+) —(SDW). , . 

[ c „ _^1.2 (<!> W ) + /„,6 ^1,1 W» W )J ( 35 ) 

In the charge-density-wave and spin-density-wave sectors, the symmetry of the dominant correlations is determined 
by finding the largest eigenvalue from the diagonalization of the 2x2 blocks 



(37) 



for r? = CDW, SDW. 

We have plotted in Figs. Q]and[5]the evolution of the response functions for arrays with N = 48 and 76, respectively, 
taking as a representative the largest eigenvalue for each kind of order, i.e. SS, TS, CDW and SDW. The dominant 
instability is found in the superconducting channel with d x 2_ y 2 symmetry. The corresponding eigenvector of the 

(SS) 

matrix R a b represents the amplitude of the Cooper pair wave function along the Fermi line, which turns out to be 
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FIG. 4: Plot of the dominant response functions for each different kind of order as a function of / = — log(A/_B c ). From 
top to bottom, we have SS d x 2_ y 2 (solid), CDW (dotted), SDW (dashed) and TS (dotted-dashed) response functions. The 
data correspond to a Hubbard model with N = 48, t' = —0.2t and U = t, at the Van Hove filling. The point where the 
renormalization group flow diverges is at I* ~ 6.96. The vertical line at l c \ ~ 6.40 indicates the energy scale at which the Fermi 
line starts to collapse. 
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FIG. 5: Similar plot as in Fig. Hfor N = 76, with I* « 4.99 and l cl w 4.36. 

peaked (with opposite signs) at the two hot spots around (jr, 0) and (0, 7r), and vanishing elsewhere. This result is 
consistent with the dominance of the d-wave correlations observed in TV-leg Hubbard ladders for small values of N. 

At this point, we can go beyond the renormalization group solution of the model, by incorporating the nonpertur- 
bative methods used in the previous section for the analysis of the quasiparticle properties. Our approach consists in 
computing the string of anomalous scaling dimensions {7a} with the effective couplings at each renormalization group 
step. To this end, we use the couplings calculated by solving the renormalization group equations to diagonalize the 
Hamiltonian in Eq. ((SJ) for < I < I* . On technical grounds, this amounts to bosonize the electron system at different 
values of A. We end up in this way with an anomalous electron dimension j a (I) that depends on the energy scale. 
The numerical computation of J a (l) reveals that, at a scale l c \ < I*, the anomalous dimension for the Fermi points 
nearest to (tt, 0) and (0,7r) diverges (namely for k a = ir/N and k a = ir — ir/N), while 7 a (Z c i) remains small for the 
rest of the electron fields. 

The divergence of the anomalous dimension j a (l) near (tt, 0) and (0, n) can be traced back to the singular behavior of 
the change of variables P3|) - ([Ht|) at the point I = l c i- The corresponding operators B a (q) and B a (q) have projections 
s a b, S-T)i t and tab that diverge for particular components B\(q) and B {q). These are made of a C4y-symmetric 
combination of operators peaked at the points nearest to (71", 0) and (0, tt). When rewritten in terms of the density 
operators B a {q) + Ba(—q), the change of variables leads however to a projection s a b + i^b (= s-r + t^) that vanishes 
at / = l c \ for the mentioned component. This shows that the situation is similar to that found in a ID system 
with two subbands, where under certain conditions the Luttinger liquid parameter K + in the channel of the total 
charge density may vanish at a certain point of the renormalization group flow22. The correspondence is such that 
Sab ~ {yjK+ + l/y/K+)/2 and tab ~ (y/ K + — l/y/K+)/2. Specifically, in the model of two subbands with equal 
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Fermi velocity v a = = vp, 



K + = 



\ 



^+(/^ ) +/g)-(/it ) +/iV) 
™ F +(f®+f$)+(fSt+f%) (38) 



We have checked that, as well as in the system of two subbands, the vanishing of s a b+tab takes place in our model when 
the interband interactions / ~t between the points nearest to (tt, 0) and (0, tt) grow large enough upon renormalization 
to drive to zero the corresponding Luttinger liquid parameter. 

We have then identified the origin of the singular behavior in the change of variables (fT5|) -(fT6 | . This takes place 
when the couplings / "t become of the order of the Fermi velocity v a , that is, when they are leaving the weak-coupling 
regime. Thus, the divergence in the anomalous scaling dimension at I — l c \ is found before entering the regime of large 
growth of the renormalized couplings, and it may be considered a more reliable result than the predictions obtained 
from the large growth of the response functions. 

The above discussion serves also to clarify the physical meaning of the divergence of 7„. This is the reflection of a 
local instability in the string of Fermi points, since at this stage it only affects the points nearest to (tt, 0) and (0, tt). 
We can first analyze the effect of the divergence of 7 Q by recalling that this anomalous dimension dictates the behavior 
of the observables defined from the electron propagator. In particular, the density of states, that we define here as a 
quantity n a (e) depending on the Fermi point, is obtained by taking the Fourier transform of the time-dependence of 
the electron propagator 

(C0M)Vw(z,*'))~ u t m +7a (39) 

Thus, we get the power-law behavior 

n„(e) ~ e 7 " (40) 

where the energy is measured in units of the high-energy cutoff. Very large values of the exponent "/ a indicate a very 
strong suppression of the electron quasiparticles, as the shape of the density of states becomes almost flat at low 
energies. We see that, at the point where j a diverges, a gap opens effectively in the spectrum of fermion excitations, 
since the density of states vanishes then over a range of low energies. 

The opening of a gap near (tt, 0) and (0, tt) in the string of Fermi points can be studied more precisely in our 
renormalization group framework. For this purpose we deal with the angle-resolved quasiparticle weight Z a , which is 
related to the anomalous dimension j a through the scaling equation 

1 dZM >--,M. (41) 



Z„(I) SI 



We have integrated numerically Eq. (|4"Tj) for every a in the range < I < l c \ . The result is shown in Fig. [5] (top 
line) for N = 48. We observe that the quasiparticle weight vanishes at I — l c \ for the Fermi points nearest to (tt, 0) 
and (0, tt), while in the rest of the Fermi line the renormalization is quite mild. We remark that the vanishing of Z a 
indicates the disappearance of fermion excitations at the Fermi point a below an energy given by Ai = E c e cl . In 
this respect, the physical effect is different than in a Luttinger liquid where the anomalous dimension of the electron 
is independent of the energy scale. A constant value of 7 implies a quasiparticle weight Z a ~ exp(— l^ a ), vanishing 
only at zero energy. In our case, however, the quasiparticle weight vanishes at an energy Ai ^ 0, implying the opening 
of a gap below that energy scale at the points near (tt, 0) and (0, tt). 

It is interesting to follow the trend of the rest of Fermi points by further lowering the energy scale. This is achieved 
by solving again the renormalization group equations for the set of Fermi points obtained by excluding the decoupled 
ones. For N — 48 we remain with 96 — 8 = 88 Fermi points, and the new initial conditions are set by evaluating 
the corresponding couplings at l c %, as obtained from the previous renormalization group flow. According to the 
above discussion, all of them are still in the weak coupling regime, making the calculation well-posed. The new 
renormalization group flow is obtained numerically starting from l c \. Again an instability is observed, due to the 
growth of the interactions involving scattering between Fermi points with smaller Fermi velocity (i.e. at k a = Ztt /N 
and k a = tt — Ztt/N). Following the same procedure shown above, the new anomalous dimensions 7 a and the 
quasiparticle weights Z a are computed. Also in this case, the j a for the Fermi points nearest to the hot spots diverge 
at l C 2 > l c \- The new quasiparticle weights are obtained by solving Eq. (|41[) in the range (Z c ii^c2), with initial 
conditions dictated by the previous integration of the same equations. In Fig. [5] (second line from top), we see that 
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FIG. 6: Angle-resolved quasiparticle weight Z a as a function of the transverse momentum k a for the first 4 iterations of the 
combined renormalization group and bosonization approaches for an array with N = 48. From top to bottom, Z a is shown at 
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FIG. 7: Plot of the critical scales I* (dotted line) and l c \ (full line) as a function of N. The parameters are the same as in Fig. 
II 

Z a vanishes again only for the points with smaller Fermi velocity. We observe that its renormalization is now more 
pronounced at different angles with respect to the previous stage. 

This procedure can be iterated, giving rise to the progressive destruction of the Fermi line, from (tt, 0) and (0, 7r) 
towards the diagonals of the Brillouin zone. It should be noted that the energy scale A„ = E c er lan at which the n-th 
Fermi point disappears may become very small. In Fig. [5] we have plotted the angular distribution of Z a for the first 
4 iterations, corresponding to the destruction of the Fermi points having transverse momentum ±(2n + l)7r/48 and 
±7r =F (2n + l)?r/48 , with n — 0, 1, 2, 3. We observe that, by increasing the number of iterations, the renormalization 
of Z a away from the hot spots becomes more and more pronounced. The plot is actually reminiscent of the results 
for the quasiparticle renormalization shown in Ref. l2Ct There is a clear indication that, while the collapse occurs on 
finite arcs around (tt,0) and (0, 7r), the rest of the Fermi line results renormalized by the interactions. 

The comparison of the plots in Figs. 2] and [5] gives a fair idea about the dependence of the above results on the 
size of the array of chains. The <i-wave superconducting correlations turn out to be always dominant (see Fig. O, 
but the collapse of the Fermi line at the hot spots prevents them from growing large in any case. In Fig. [7] we have 
represented the trend of the critical scales I* and l c \ as a function of N. It appears that I* approaches an asymptotic 
value in the large- N limit, while the difference between I* and l c \ does not vanish by increasing the size of the array. 
This suggests that the development of the pseudogap should survive in the 2D limit, ending up with a paradigm of 
electron liquid in which the Fermi surface is progressively destroyed as the system is probed at lower energy scales. 



IV. CONCLUSION 



In this paper we have studied the low-energy electronic instabilities in a 2D Hubbard model formed by a large array 
of coupled chains. Within a low-energy cutoff E c about the Fermi level, the different interactions can be classified into 
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a reduced number of channels^ -, which coincide essentially with those obtained from the analysis of the kinematics 
for a 2D Fermi surface^ . In the model with bare on-site repulsion U, the different couplings in momentum space have 
a strength of the order of U /N, which makes appropriate the use of renormalization group methods to determine the 
scaling of the couplings as the cutoff is reduced towards the Fermi level. 

Complementing this approach, we have also applied bosonization methods to describe the dressing of the electron 
fields at the different Fermi points. Thus, we have computed the string of anomalous scaling dimensions {7a} with 
the effective couplings at each renormalization group step. This combines the power of the renormalization group 
and the bosonization methods, allowing to characterize the points where the metallic phase may break down. In this 
way, we have been able to discern the competition between the tendency to the formation of ordered phases in the 
2D electron liquid and the charge instabilities that modify the shape of the Fermi line. 

Our main interest has been focused on the instabilities that appear when the Fermi level is close to the Van 
Hove singularity arising from the saddle points at (7r,0) and (0,7r). It is well-known that the singularity in the 
density of states leads to strong correlations that signal the tendency towards a spin-density-wave phase or d-wave 
superconductivity at low energies. This has been studied in the past using mainly renormalization group methods. An 
important drawback in this approach comes however from the unconventional form of the 2D perturbative expansion, 
with terms that diverge at low energies as the square of log(e) and which are not therefore properly renormalizable. 
This may cast some doubts about the complete predictability of the renormalization group to discern the competition 
between low-energy instabilities in the 2D model. 

Alternatively, it has been pointed out that the renormalization of the quasiparticle properties has to be quite 
significant near the Van Hove singularity, and that this may lead to a drastic softening of the electron correlations^ -. 
In our approach, we have recovered some of the trends found within the conventional renormalization group schemes, 
like the strong renormalization of the dispersion relation around the saddle point o 15 i 23 or the growth of the d-wave 
superconducting correlation a 11 ! 12 ] 14 ! 16 ] 17 ! 18 ] 19 ! 21 ! 22 . However, our main finding is that, at the Van Hove filling, the 
progressive destruction of the Fermi line takes place around the saddle points before any correlation may diverge in 
the low-energy limit. This is a sensible result in our model, as the ID renormalization group approach is well-defined, 
while the analysis of the breakdown of the metallic phase is carried out with a nonperturbative bosonization approach. 

This study may have also implications for the understanding of the pseudogap phase in the copper-oxide super- 
conductors. The loss of fermion excitations takes place in our model as the electron fields are attenuated due to the 
dressing from the renormalized interactions at the hot spots around (n, 0) and (0,7r). This happens as the charge 
stiffness goes to zero at some of the Fermi points in the low-energy limit. We stress that this behavior marks the 
approach of a critical point, characterized by the vanishing of some of the Luttinger liquid parameters encoded in the 
expression (|19p for the anomalous dimensio n 32 ! 33 . The critical point is reached as some of the effective couplings leave 
the weak-coupling regime, so our description of the pseudogap does not rely on strong-coupling features to account 
for the destruction of the Fermi line. In future studies, it would be interesting to investigate the extrapolation of 
the present results to systems with larger number of chains. It seems plausible that the destruction of the Fermi line 
we have described should persist for increasing number of chains, as it relies on a critical behavior that is a robust 
property of the bosonization of the system in the low-energy limit. 
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